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ABSTRACT 

ILoeb k. Zaldarriagal ((2003) have recently proposed that observations of the CMBR 
brightness temperature fluctuations produced by the HI inhomogeneities prior to reion- 
ization hold the promise of probing the primordial power spectrum to a hitherto un- 
precedented level of accuracy. This requires a precise quantification of the relation 
between density perturbations and brightness temperature fluctuations. Brightness 
temperature fluctuations arise from two sources (1.) fluctuations in the spin tempera- 
ture, and (2.) fluctuations in the HI optical depth, both of which are caused by density 
perturbations. For the spin temperature, we investigate in detail its evolution in the 
presence of HI fluctuations. For the optical depth, we find that it is affected by density 
perturbations both directly and through peculiar velocities which move the absorption 
features around in frequency. The latter effect, which has not been included in earlier 
studies, is similar to the redshift space distortion seen in galaxy surveys and this can 
cause changes of 50% or more in the birghtness temperature fluctuations. 

Key words: cosmology: theory - cosmology: large scale structure of universe - diffuse 
radiation 



1 INTRODUCTION 

The possibility of probing the universe at high redshifts 
using the HI 21 cm line has been the topic of exten- 
sive theoretical investigation. This is perceived to be the 
most promising window for studying the "dark ages", 
the era between the decoupling of the CMBR from the 
primeval plasma at z ~ 1000 and the formation of 
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proaches T 7 . This gives a range of redshifts where T 3 < T 7 . 
We then have a window in redshift 30 < z < 200, or equiva- 
lently in frequency v = 1420MHz/ (1 + z) where the HI will 
produce absorption features in the CMBR spectrum. 



the first luminous objects at 
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ture T g is maintained at CMBR temperature T 7 through 
collisions of the CMBR photons with the small fraction of 
electrons that survive the process of recombination. The col- 
lision process becomes ineffective in coupling T g to T 7 at 
2 ~ 200. In the absence of external heating at z < 200 the 
gas cools adiabatically with T g oc (1 + z) 2 while T 7 oc (1 + z). 
The spin temperature T s is strongly coupled to T 9 through 
the collisional spin flipping process until z ~ 70. The colli- 
sional process is weak at lower redshifts, and T s again ap- 
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In a recent paper ILoeb fc Zaldarriagal ([2003) propose 
that observations of the angular fluctuations in the CMBR 
brightness temperature T;, arising from the HI absorption 
can be used to study the power spectrum of density fluctua- 
tions at small scales to a level of accuracy far exceeding those 
achievable by any other means. The enormous wealth of in- 
formation arises from the fact that observations at differ- 
ent frequencies which are sufficiently separated will provide 
independent estimates of the power spectrum at the same 
wave number k. These observations will probe the power 
spectrum before the epoch of structure formation, and they 
hold the possibility of revealing the entire primordial power 
spectrum down to very s mall scales. In another recent paper 
ICnedin fc Sheveil l)2003h have studied the linear fluctuations 
in the 21 cm emission from the pre-reionization era. They 
show that that it should be possible to detect this signal 
against the foreground contaminations in the frequency do- 
main. This signal is expected to constrain the equation of 
state of the universe at high z. 
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The CMBR brightness temperature is related to T s 
and the HI number density tih as Tj oc (1 — T 7 /T s ) uh- 
Fluctuations in T& arise from fluctuations in rm directly 
and also through fluctuations in T 3 which in turn arise 
fro m fluctuations in nn - In ca lculating the fluctuations in 
T 3 , lLoeb fc Zaldarriaeal (l2003Tl consider only one process, 
namely the change in the collision rate arising from fluctua- 
tions in riH. Perturbations in riH will also produce perturba- 
tion in T g , which in turn will affect T s . This effect has not 
been taken into account in their work. 

Density perturbations produce peculiar velocities. This 
causes the frequency of the HI absorption features to be 
shifted by the line of sight component of the peculiar ve- 
locity. This effect will rearrange the HI absorption features 
in frequency space where converging velocity patterns ap- 
pear as enhancements in the HI number density and di- 
verging velocity patterns appear as decrements in the HI 
number density. It may be noted that this i s the fa miliar 
linear redshift space distortion (Kaiser effect, lKaiseiiri987h 
seen i n galaxy redshift surveys. This effe ct has been stud- 
ied bv lBharadwai. Nath and Sethil (l200lT) in the context of 
cosmological HI emission from z ~ 3.5. 

It is important to identify and take into account all 
possible contribution to the brightness temperature fluctu- 
ations, if these are to be used to extract precise information 
about the power spectrum and the equation of state of the 
universe. In this paper we study two effects which will con- 
tribute to brightness temperature fluctuations, namely (1.) 
perturbations in the gas temperature produced by density 
fluctuations, and (2.) the effect of redshift space distortions. 
To the best of our knowledge, these effects have not been 
included in earlier works. 

We next present an outline of the paper. In Section 
2. we discuss the processes involved in determining the 
brightness temperature fluctuations and present the rele- 
vant equations. In Section 3, we present our results and dis- 
cuss their consequences. It may be noted that we have used 
(n m o, Oao, Qboh 2 , h) — (0.3,0.7,0.02,0.7) whenever specific 
values have been needed for the cosmological parameters. 



2 CALCULATING THE BRIGHTNESS 
TEMPERATURE FLUCTUATIONS. 

We first consider the evolution of the gas temperature after 
the recombination era (z ~ 1000) when it becomes largely 
neutral. This is governed by the equation 



(1 + Z ) 3/2 (T 7 -T 9 ) (1) 



dT g _ 2Tg_dnu_ _ -9.88 x 10~ 8 
dz 3riH dz Q.^h 2 

The third term in the above equation represents the energy 
transfer from the CMBR t o the gas thr ough the collisions 
with the residual electrons <IPeebleall993r) . This terms tries 
to maintain the gas temperature at the CMBR temperature 
as the universe expands. The second term is the change in 
T g in adiabatic expansion. If the HI is uniformly distributed, 
then nH oc (1 + z) 3 and in the absence of the CMBR heating 
we have T g oc (1 + z) 2 . Collisions are able to maintain T g = 
T 7 = 2.73 K(l + z) up to a redshift z ~ 200 (Figure[TJ after 
which T g oc (1 + zf. 

We next consider the evolution of the spin temperature 
T s which is defined through the relation 
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Figure 1. This shows the evolution of the CMBR temperature, 
the gas temperature and the spin temperature as the universe 
expands. 
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where m, no are the population densities and g\ — 3, go — 1 
the spin degeneracy factors of the excited and the ground 
states of the HI 21 cm transition. In this equation — 
hpVe/kB = 0.068K, where hp is the Planck's constant, v e = 
1420MHz is the frequency corresponding to the 21 cm line 
and fcfl is Boltzmann's constant. The evolution of the ground 
state population density is governed by two processes, one 
collisional and the other radiative 



(l + 4) 



no = mCio — noCoi 
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where a(t) is the scale factor, Coi and Cio are the collisional 
excitations and de-excitation rates of the hyperfine levels, 
Aio is the Einstein spontaneous emission coefficient, Boi 
and Bio are the Einstein B coefficients and 7„ e is the specific 
intensity of the background radiation at v £ . 

In the regime of our interest T B ,T g ,T n 3> T* and we 
can use the approximation e~ T * //T = 1 — (T+/T) through- 
out. Also, the fact that in equilibrium the collisional pro- 
cesses and the radiative process are separately balanced 
gives us the relations Coi = 3(1 — T- k /T g )Cio and Boi = 
3B 10 = (3\ 3 rJ2h v c)A w w here A w = 2.85 x Hr 15 s -1 
jRvbicki fc Lightmanl fl979h . The collisional de-excitation 
rate can be written as Cio = (4/3)«(l — 0)nH where 
the values of k(\ — 0) i s tabulated as a function of T g 
( Allis on fc Dalgarnolll969l) . Using these and equation (|3J we 
obtain an equation for the redshift evolution of T s 



— ( — ) 

dz \T S J 



T. 



H(z){l + z) 
1 

~ Ts 



Cio + 



-Ai 



(4) 



where H(z) is the Hubble parameter. The collisional term 
tries to set the spin temperature at the same value as the gas 
temperature while the radiative term tries to set it at the 
CMBR temperature, which process dominates being decided 
by the rate coefficients. At high redshifts the collisional pro- 
cess dominates and the spin temperature closely follows the 
gas temperature. At lower redshifts «h falls substantially, 
the collisional process looses out to the radiative process 



The CMBR fluctuations from HI perturbations prior to reionization 3 



0.6 
0.5 
0.4 
0.3 
0.2 
X 0.1 
M o 
-0.1 
-0.2 
-0.3 



Figure 

defined 



g(z) 
s(z) 




100 

z 



2. This shows the evolution of the functions g(z) and s(z) 
in the text. 



and the spin temperature approaches the CMBR tempera- 
ture. Figure shows the evolution of the spin temperature 
as the universe expands. 

We next shift our attention to the effect of HI density 
perturbations Ah(x,z) = (iih(x,z) — riH(z))/riH(z). These 
will produce fluctuations in the gas temperature. If the gas 
were undergoing adiabatic expansion, the fluctuations in the 
gas temperature A s (x, z) = (T g (x, z) — T g (z))/T g (z) would 
be related to Ah through A s = (2/3)Ah. This will get mod- 
ified because of the energy that is pumped into the gas from 
the CMBR and A 9 = during the era when T g = T 7 . We 
define a function g(z) = such that A g (z) = g Ah(z). 

Using this in equation Q we obtain 



dg 9.88 x 10" 8 T 7 



dz 



n b h 2 T a 



3/2 



1 dAn 
A H dz 



(5) 



The first term in RHS arises from the coupling of the gas 
to the CMBR and it tries to set g[z) — > while the second 
term corresponds to adiabatic expansion and it tries to make 
g(z) — > 2/3. The quantity g(z) is expected to evolve from 
g(z) — to g(z) — 2/3 as the universe expands and the 
contribution of the heat pumped into the gas goes down. 

An interesting feature is that g(z) depends on the 
growth rate of density fluctuations. For example, it follows 
from equation that a static density perturbation which 
does not evolve in time will not produce fluctuations in 
the gas temperature. Here we assume that Ah follows the 
dark matter perturbation and grows as Ah oc a(z). The 
result of integrating equation © is shown in figure [5] We 
see that g(z) ~ 0.3 in the redshift range of our interest 
(30 < z < 100). 

We finally come to the fluctuations in the spin tem- 
perature A 3 (x,z) = (T s (x,z) — T s (z))/T s (z) produced 
by density perturbations. From equation @ we see that 
fluctuations in T 3 can arise from changes in T g and from 
changes in the collision rate. The changes in collision rate 
Cio = (3/4)k(1 — 0)riH will come about directly through 
changes in nu and also through changes in T g which will 
affect the value of k(1 — 0). Taking into account both these 
effect we have ACio = (l + (2/3) dlnn/dlnT g ) Cio A H . Defin- 
ing a function s(z) = §^-. such that A 3 (z) = sAh(z) and 



using equation Q we obtain 
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Here again, the evolution of s(z), like that of g(z), de- 
pends on the time evolution the density fluctuations. Fig- 
ure |21 shows the evolution of s(z) under the assumption 
Ah oc a(z). We find that s(z) > for z > 90 ie. a posi- 
tive density perturbation causes the spin temperature to go 
up, and the effect is opposite at z < 90. 

During the era when T s < T 7 the HI along a line of 
sight n reduces the CMBR brightness temperature at the 
frequency v by an amount 



T b (n,u) = (T s - T T ) r/(l + z) 



(7) 



Here r is the optical depth of the 21 cm HI absorption given 
by 



3 nu hp c 2 Aio a 2 (z) . dr 
32 7rfcg7> e ' 



(8) 



where r is the comoving distance to the HI whose 21 cm 
absorption is redshifted to v. 

We are interested in the angular fluctuations of the 
brightness temperature Tt,(n, v). HI density fluctuations will 
produce fluctuations in Tf,(n, v) through the fluctuations in 
the spin temperature discussed earlier. Density fluctuations 
will also directly affect 21, through variations in the opti- 
cal depth which we now calculate. The relation between the 
comoving distance r and the frequency v is given by 



cda 
a 2 H(a) 



(9) 



where v is the line of sight component of the peculiar velocity 
of the HI which produces the absorption. Density perturba- 
tions will, in general, be accompanied by velocity perturba- 
tions and these will move around the HI absorption features 
in frequency. Here we assume that the HI traces the dark 
matter and that we can use linear perturbation theory to 
relate the peculiar velocity to the density perturbations. In- 
corporating the effect of both the density fluctuations and 
the peculiar velocity we have 



3tih hp c Aw 
32 it kpTsvi H (z) 



Au- 



di' 



H{z)a{z) dr 



(10) 



Here we have dropped terms of order v/c in the final expres- 
sion. Also, we have retained terms only to linear order in v. 
There is also the effect of our own motion which we have 
not included. These effects are not expected to be impor- 
tant. The term involving the derivative of the peculiar ve- 
locity is the dominant effect, particularly at the small scales 
of interest here. 

Combining the effects of the fluctuations in the optical 
depth and in the spin temperature we can write the fluctu- 
ations in the CMBR brightness temperature as 



ST b (n,v)=T 



where 



(a h 



1 dv 
H a dr 



-sAt 



(11) 
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It is convenient to express this in Fourier space where 

AhM= / (03 e_<k ' XA ( k ^) ( 13 ) 

and the Fourier transform of the peculiar velocity is given 
by v(k,z) = -iif (z)o(z)kA(k, z)/k 2 . Using this we express 
the fluctuations in brightness temperature as 

d 3 k 



<m(r 



= T 



-ikrfi 



(2tt)3 
T 3 



A(k,z) 



(14) 



where [i is the cosine of the angle between the comoving 
wave vector k and the line of sight n. 

We now calculate the angular power spectrum of 
the brightness temperature fluctuations resulting from the 
density fluctuations A(k, z). The statistical properties of 
A(k, z) are specified through the 3D power spectrum de- 
fined as 



<A(k,z)A(k ,z)) = (27r) 3 (£(k-k )P(k,z) 



(15) 



where (...) denotes ensemble average and Sd() is the Dirac 
delta function. 

The angular power spectrum is calculated by decompos- 
ing the angular dependence of 8Tb into spherical harmonics 
with expansion coefficients ai m (y) and using these to cal- 



culate the angular power spectrum Ci(y) 



(I In 



2 ). The 



angular power spectrum can be expressed in terms of the 
3D power spectrum as 



C ! (i/)=4ttT 2 



(2^)3 



P(k,z) 



Ji(kr) 



+ Y s 3i( kr ) 



(16) 



where ji(kr) are the spherical Bessel functions and 
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3 RESULTS AND DISCUSSION. 

In this paper we have investigated in detail the CMBR fluc- 
tuations produced by HI priot to the epoch of reionization. 
As proposed by Loeb and Zaldarriaga (2003), this holds the 
promise of allowing the power-spectrum of density fluctua- 
tions to be probed to a high level of precision. 

HI density perturbations produce fluctuations in the 
decrement of the CMBR brightness temperature by two 
means (1.) through fluctuations in the optical depth, and 
(2.) through fluctuations in the spin temperature. The ef- 
fect of changes in the optical depth in response to a pos- 
itive density perturbation (curve A of Fig. [3J is such that 
it reduces 21, and enhances the decrement in the brightness 
temperature. This effect is maximum at z ~ 80. This effect 
is enhanced by peculiar velocities. 

The change in brightness temperature produced by HI 
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Figure 3. This shows how the CMBR brightness temperature 
fluctuations (in units of T) respond to HI density perturbations. 
A. shows (1 — Tj/Ts) which is how it responds through changes 
in the optical depth, ignoring the effect of peculiar velocities. Pe- 
culiar velocities will enhance this effect. B. shows sT 7 /T s which 
is how it responds through changes in the spin temperature. C. 
shows the total response for pu 2 = 2/3. D is the same as B without 
the gas temperature fluctuations. 



density perturbations through changes in the spin tempera- 
ture varies with z (curve B of Fig. El - Density perturbations 
increase the spin temperature and the brightness temper- 
ature in the redshift range z > 100. Here the collisional 
process is very efficient and Ts closely follows T g . A positive 
density perturbation increases T g which causes T s to also in- 
crease. The effect of density perturbations on the brightness 
temperature acting through changes in the optical depth and 
through the spin temperature have opposite signs. The effect 
of changes in the optical depth is larger and the brightness 
temperature decrement is enhanced by a positive density 
perturbation (curve C of Fig. - 

In the redshift range z < 100 positive density per- 
turbations lower the spin temperature which enhance the 
brightness temperature decrement. In this regime the colli- 
sional process slowly looses out to the radiative process and 
T s — > T 7 . A positive density perturbation enhances the colli- 
sion rate which pulls the spin temperature down toward the 
gas temperature. The two processes which contribute toward 
brightness temperature fluctuations both act to enhance the 
decrement. Curve C of Fig. [3] shows the combined effect of 
both these processes for the value /i 2 = 2/3. We find that 
the response of the brightness temperature to density per- 
turbations peaks at z ~ 55. This is somewhat sm aller than 
the value obtained bv lLoeb fc Zaldarriaeal ([2003). Curve D 
of Figure 3 shows the contribution to brightness tempera- 
ture fluctuations arising from the spin temperature if the 
effect of density perturbatio ns on the gas temperatu re is 
not taken into account (eg. lLoeb fc Zaldarr iaga 2003) . We 
find that including the gas temperature makes a significant 
change, particularly at z > 100 where there is a qualitative 
difference between curves B and D. 

We have calculated the angular power spectrum Ci (v) of 
the brightness temperat ure fluctuations for th e COBE nor- 
malized ACDM model llBunn fc Whitdll996ft . The power 
spectrum has been suppressed beyond the arbitrarily cho- 
sen value k = 14/iMpc -1 using a Gau ssian cut-off. Our re- 
sults a re in qualitative agreement with lLoeb fc ZaldarriagJ 
( 2003). We find that the signal peaks at z ~ 50 (Fig. |IJ 
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Figure 4. This shows the angular power spectrum of the CMBR 
brightness fluctuations at various redshifts with and without the 
effects of peculiar velocities. 



where the product of the growing mode of density perturba- 
tions and the response of brightness temperature to density 
perturbations is maximum. 

To gauge the effect of peculiar velocities, we have calcu- 
lated Ci(y) without incorporating this effect. This is easily 
done by replacing Ji(kr) with ji(kr) in equation (I16t . The 
results are shown in Fig. [I] We find that the peculiar veloc- 
ities increase y/Ci by more than 50%. 

We have also quantified the effect of gas temperature 
fluctuations. We find that for z < 100 the values of \fC[ 
are ~ 10% lower if the gas temperature is not taken into 
account, and the effect is reversed at z > 100 where \fCi is 
~ 30% higher. 

The ability to probe the dark matter power spectrum 
using the Ci(v)s will be restricted to scales k < fcj where 
fcj is the Jeans wave number. This has a nearly constant 
value ~ 500/iMpc -1 in the redshift range of interest. The 
HI power spectrum on scales smaller than the Jeans length 
scale is interesting in its own right. The HI perturbations will 
undergo acoustic oscillations on these scales. The spin tem- 
perature fluctuations and the peculiar velocities produced by 
density perturbations will be quite different from the situa- 
tion considered here. On scales slightly larger than the Jeans 
lengthscale (2n/kj) the density nuctuations are mildly non- 
linear with rms. values in the range .7 > a > .1, and the 
Z'eldovich approxim ation may give a better description (eg. 
iHui fc Gnedrnlll997l) . 

The low frequency cut-off imposed by the earth's iono- 
sphere restricts observations to frequencies more than ~ 
10 — 20 MHz. Extracting the HI signal from the contam- 
inations arising from the galactic and extragalactic fore- 
grounds is going to be a big challenge. The foregrounds are 
mostly continuum sources whose contribution varies slowly 
with frequency. It will be necessary to combine both the 
angular fluctuations and the frequency domain properties 
of the CMBR br ightness temperature fluctuations in order 



of the IMHK brightness temperature fluctuations m order 
to detect it fee. IShaver et al.lll99l IDi Matteo et alJl200l 
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